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^ ' Abstract 



The tunneling effect of a periodic potential with an asymmetric twin barrier 
per period is calculated using the instanton method. The model is derived 
from the Hamiltonian of a small ferromagnetic particle in an external mag- 
netic field using the spin-coherent-state path integral. The instantons in 
two neighbouring barriers differ and lead to different level shifts Ae 1; Ae2- 
We derive with Bloch theory the energy spectrum which has formally the 
structure of an energy band. The spectrum depends on both level shifts. 
The removal of Kramer's degeneracy by an external magnetic field is dis- 
cussed. In addition we find a new kind of quenching of macroscopic quantum 
coherence which is irrelevant to Kramer's degeneracy. 

PACS numbers: 11.15.KC, 03.65.Db, 11. 10. St, 73.40.GB 
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1. Introduction 



Quantum tunneling with various solvable potentials such as the double-well 
Ufl, inverted double-well 0] and sine-Gordon || potentials have been stud- 
ied extensively in the framework of the instanton method. Instantons in field 
theory of 1 + dimensions are viewed as pseudoparticles with trajectories 
existing in barriers, and are therefore responsible for tunneling. Since in- 
stanton trajectories have zero energy and nontrivial topological charge, and 
interpolate between degenerate vacua, the instanton method is only suitable 
for the calculation of tunneling between neighbouring vacua. Recently a new 
type of pseudoparticle configuration has been constructed ||, |j and dubbed 
"periodic instanton". These periodic instantons have finite, nonzero energy 
and are found to be useful in the calculation of tunneling at excited states 

In the case of the double-well potential the amplitude for transitions be- 
tween degenerate vacua induced by instantons results in the level shift of the 
degenerate ground state. In the case of the periodic potential the instanton 
induced transition amplitude is valid for tunneling through only one of the 
barriers and is the same for all barriers. The tunneling amplitude results 
in the level splitting which extends to an energy band || |6| by translation 
symmetry, as follows from Bloch theory. 

The present paper is devoted to solving the tunneling problem with a 
periodically recurring asymmetric twin barrier which leads to two kinds of 
instantons. This case is of interest in its own and worthy of development of 
the appropriate instanton technique, and is also of practical importance for 
macroscopic quantum tunneling J7| in magnetism j8|. 

The layout of the paper is as follows. In sections 2 and 3 the effective 
Lagrangian is derived from the Hamiltonian of a small ferromagnetic par- 
ticle in a constant magnetic field. The exact instanton solutions existing 
in the barriers are given. We also calculate the total action along the in- 
stanton trajectories (the socalled instanton mass). Section 4 is devoted to 
the calculation of the energy spectrum with the help of Bloch theory. The 
resulting formula involves two level shifts to be determined from tunneling 
amplitudes. In section 5 the tunneling amplitudes are obtained with the 
instanton method up to and including the one loop correction. Finally we 
discuss the quenching of macroscopic quantum coherence (MQC), i.e. that of 
the energy band due to the topological phases of the effective wave functions 
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of the small ferromagnetic particle in the magnetic field. 



2. The effective Lagrangian of a small ferromagnetic 
particle in a constant magnetic field 



We consider a giant spin in a constant magnetic field described by the 
Hamilton operator 

H = Kis 2 g + K 2 s 2 y - 2fi B Bs y (1) 

where s^i = x,y,z, are spin operators satisfying the usual commutation 
relation [si, Sj] = ie^Sk (using natural units throughout), and hb is the 
Bohr magneton. The model describes in the absence of the magnetic field B 
XOY easy plane anisotropy |9[] and an easy axis along the x-direction with 
the parameters K\ > K 2 > 0. The external magnetic field is assumed to 
be perpendicular to the easy direction. Our purpose in the following is to 
convert the quantum spin system into a potential problem suitable for the 
study of quantum tunneling. To this end we start from the coherent state 
representation of the matrix element of the time-evolution operator with the 
Hamiltonian of eq. ([!]). With the help of the spin-coherent path integral we 
obtain 

< n/ |e- 2 ^ T |n, >= e-^-^'K^tf^M (2) 

where 

Kfo, t f ; fa U) = J V{^V{v\e^ C{Mdt (3) 

is the path integral in phase space with canonical variables 4> and p = s cos 9. 
Also 

C = 4>p-H{4>,p) (4) 
is the phase space (or first order) Lagrangian. The Hamiltonian 

has position dependent mass m(<f>) and potential 

m ^ ~ 2#i(l - A sin 2 + aA sin 0) ' 

V{4>) = K 2 s 2 (sin 2 <p -2a sin <p + a 2 ) (6) 
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respectively. The potential minima are chosen to be zero (see Fig.l). In these 
expressions 

K 2 HbB 
K\ sK 2 

are dimensionless parameters whose values are correspondingly indicative of 
the degree of anisotropy of the ferromagnetic particle and the asymmetry of 
the twin-barrier (Fig.l). The kets > and |nj > denote the initial and 
final spin-coherent- states and IT = tf — ti is the time difference. Here 
s = s(sin#cos0,sin#sin0, cos 9) is visualised as a classical spin vector with 
spin number s, polar angle 9 and azimuthal angle (p. In the above derivation, 
the large spin limit s » 1 has been used since giant spins with spin quantum 
number s >> 1 are believed to correctly describe ferromagnetic grains. To 
simplify the model we also consider the tunneling near the easy plane, i.e. 
for 9 ~ ~. Our rule of simplification is then that for s — > oo,cos# — > 
the product s cos 9 is to remain finite. The phase factor e~^ f ~^ s in the 
transition amplitude eq. (@) is the characteristic of the spin system and 
can be put into the Lagrangian, i.e. the expression <fif — 0j = f t * <fidt, and 
can be identified as a Wess-Zumino term || . The effect of the phase factor 
has been studied extensively || [10| . The position dependent kinetic term 
requires that one should start from the phase space path integral [H|, i.e. 
eq.(Q). Integrating out the momentum in the path integral eq. (|3p, we obtain 
the usual Feynman propagator in configuration space, i.e. 

m,tf,M = fvifte'tt^* (7) 
with the second order Lagrangian 

C = \m^)ft - V(<P) (8) 
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3. Exact instanton solutions, instanton "mass" and re- 
duced tunneling model 



The instanton solution in 1 + dimensions is a trajectory of the pseu- 
doparticle which exists in the classically forbidden region, namely the poten- 
tial barriers, and can be obtained from the Euclidean version of the action 
of eq. (0) by Wick rotation r = it, (3 = iT. The Euclidean propagator is 



K E = J P{0}e" 5 - (9) 
with the Euclidean action defined by 

S E = jP C E dr (10) 
where the Euclidean Lagrangian is 

C E = l -mm 2 + V(<f>) (11) 

^From now on = ^ denotes the imaginary time derivative. Since there are 
two kinds of barriers we expect two types of instanton solutions. The two 
classical solutions 0w(r),i = 1,2, which minimize the Euclidean action, eq. 
(|10D, are found to be instantons with nonvanishing topological charge, i.e. 

2,/a 



u)o(ct- — 2\f\)^faZ 

.{(a 3 2 - « x 2 )n (c&W^t)), « 3 2 , k) + ai 2 F (*<%<V)), k) ) (12) 



and 



2,/a 



u {a + - 2yf\) 

2 



^ ac - « 2 2 )H ($^(#(r)), a A \ k) + a 2 2 F (*«(^(r)), k) j (13) 

where F denotes the incomplete elliptic intregral of the first kind with elliptic 
modulus k, and II is Legendre's incomplete elliptic integral of the third kind 
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with the same elliptic modulus. The corresponding parameters are defined 
by 



_ 8JA(4 + A« 2 ) 



a + a_ 



(14) 



where 
and 



4y/A 

a + 



4^A 



, ,.v- , 2(q_ - 2yf\) , 2(a + - 2yf\) 
, a 2 = — — , a 3 = — — — , a A = ; — — (15) 



a_ 



(1 — a)a + 



(1 + a)a_ 



a+ = V4 + Aa 2 + 2V\ + av^, a_ = V4 + Aa 2 + 2V\ - aV\ 
Here <&W,i = 1,2, are functions of the classical trajectories 0^(r), 



$(!) 



$( 2 ) 



sin 



sin 



-l 



(1 - sin 0^)0+ 
^2[a+-2 v / A(l + sin^ 1) )] 

(1 + sin 0^)a_ 



(16) 



(17) 



(18) 



^ 2[a_ -2^(1 -sin<^ 2) )] 

The trajectories of instantons (p^ and are shown in Fig. 1. At initial 
time Tj = — /?i (or /3 2 ) the instanton (or (/>( 2 )) starts from the minimum 
of the potential well, for instance from <pi = sin -1 a (or n — sin -1 a), and 
reaches the neighbouring minimum (pf = n — sin -1 a (or In — sin -1 a) at final 
time Tf = (3i (or (3 2 ), where 

2v^T 



cu 



(a_ - 2^)^ 



(«3 2 -ai 2 )n($«(a,A),a3 2 ,£;) 
+a! 2 F($ (1) (a,A),A;) 



$ (1) (a,A) = sin- 1 — 
« 3 



(19) 



and 



cj (a+ - 2v / A) A /a7 



(a 4 2 -a 2 2 )n($ (2) (a,A),a 4 2 ,A;) 
W^ (2) («,A),A;) 
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$ (2) (a,A) = sin- 1 — 

«4 



(20) 



It is interesting to observe that because of the position dependence of the 
mass the instantons here do not need an infinite length of time to tunnel 
through the barriers as in the usual instanton case of the double well and sine- 
Gordon potentials. The minimum action, called instanton "mass" , evaluated 
along the instanton trajectories is 



ft 



V) 2 dr = As 



A^X 



n($ (1) («,A),ai",fc) 



y/a + a 
2V\(l + a 



F($«(a,fc),fc), 



(21) 



f^fdr = As 



2^/\(a + + a_) — a+a_ 



n($ (2) (a,A),a! 2 ,A:) 



'a + a_ 
2\f\)a + + 2\f\aa. 



(22) 



For the sake of simplicity we consider in the following the extremal model 
with A << 1. This condition is of practical interest since the macroscopic 
tunneling effect is negligible for the ferromagnetic particle with large spin s 
except when A is extremely small (JEJ. The mass defined in eq.(^) may then 
be replaced by a constant, i.e. m « With this replacement the effective 
Lagrangian eq. ([TTj) reduces to a simple model. The classical solutions which 
minimise the Euclidean action are seen to be also instantons, i.e. 



2 tan 



2 tan' 



-i 



a 



i 1 



where 



/ (1) («) 



a 



a 



e VT^Mr-r ) _ /(i) ( a ) 



r Vl-a 2 ^o(r-r ) _|_ /( 2 )( a ) 
3 -v / L~^; (t-to) _ /(2)( a ) 



1- a- vr 



a 



2' 



(2), 



a 



1 + a + Vl - a 2 
1 + a — a/1 — « 2 



(23) 
(24) 

(25) 
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These solutions of the simplified model also coincide with the A — > (k — > 0) 



limit of eqs. ([T2] , p^|) . 

The trajectories of instantons 4>^\i = 1,2, are the same as those shown 
in Fig. 1, but with f3\ and (3 2 tending to infinity. r is the "position" of the 
instanton. When r = r the instantons arrive at the centres of their barriers 
((p = |, mod 2ir, for anc l = ^ 1 mod2-yr, for 4>^). The minimum actions 
(masses) evaluated along the trajectories of the instantons are seen to be 

/oo 
m^fdr = sV\{2Vl - a 2 - an + 2a sin" 1 a), (26) 
-OO 

/OO 
m(0f) 2 (ir = sv / A(2v / T^2 + a7r + 2a gin" 1 a) (27) 
-oo 

The amplitudes for tunneling through the barriers of type "1" (those of small 
size with centres at (f> = |, mod 2n) and of type "2" (those of large size with 
centres at (f) — , mod 27r) can be naively calculated by evaluating the 
Feynman propagator eq. (2.7). However, it is not clear what the physical 
result relating to the two tunneling propagators would be. In the double 
well potential case the Feynman propagator is related to the level splitting of 
the degenerate ground states. In our case we would have two different level 
splittings Aei and Ae2- To solve this apparent paradox we calculate in the 
following section the energy spectrum with the help of Bloch theory for the 
periodic potential. 



4. The energy spectrum 



We let | eo, $2n > be the ground state of the harmonic oscillator approxi- 
mated Hamiltonian i^ 1 '* in the 2nth well around the local potential minimum 
at $2n = 2rnr + sin -1 a such that 

Hj ) 1) \e ,^>=e \e ,^> (28) 

where the Hamiltonian is seen to be 

^ 1) = ^ + ^W(0-*S) 2 (29) 
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and <fi — $2n 1S the angle deviation from $ 2 ^. Similarly |eo, 3>2n+i > denotes 

* (2) 

the ground state of the harmonic oscillator approximated Hamiltonian H Q 



in the (2n + l)th well around the local potential minimum $ 2 n+i 
— sin -1 a. i^ 2 " 1 is seen to be the same as Hq but with the angle i 

(2) 

understood as the deviation from the equilibrium position $ 2n+1 . According 
to Bloch theory the true ground state of the system is a superposition of 
|e ,$ (i) >,* = 1,2, i.e. 



{2n + 



H >= E (e«*-|eb,*£ > +e^ko,$S + i >) (30) 
which is a Bloch state. £ denotes the Bloch wave vector and should be an 



integer in our case because of the 2ti boundary condition |T3[ . We then have 

H\ij>=E\^> (31) 

H is the Hamiltonian of the system with periodic potential and E is the true 
ground state energy to be determined. Simple calculation leads to (taking 
into account only nearest neighbours and remembering that as in eq.(|2|) the 
spin coherent state representation leads to wave functions with topological 

phase, e. g. (0|e o , $2} = ~ CV ( ^" )s ) 

E = e + Aei cos (s + £)(7r — 2 sin -1 a) + Ae 2 cos (s + £)(7r + 2 sin -1 a 
where 



Ae 2 



(2) 



(32) 
(33) 



(2) 



(34) 



are just the level shifts due to tunneling through the small (type 1) and 
larger (type 2) barriers respectively and they are the same for all n. The 
function %(</> — $W) denotes the ground state wave function of Hamiltonian 
F (i) . The spin number s in eq. ( |3"2"| ) comes from the inner product of the 
spin-coherent-states like the phase factor in eq.(^). The expressions Aei and 
Ae 2 can be evaluated with standard instanton methods, as we shall show in 
the following. 
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5. Amplitude for tunneling through the asymmetric 
twin barrier 



We start from the Euclidean version (imaginary time) of the element of 
the evolution operator, eq.(H), that is, 

< rvle-^ln, >= e^^'/C^,, r f ; n) (35) 

Inserting the complete set of eigenstates of the harmonic oscillator approxi- 
mated Hamiltonians and on the left hand side of eq. (|35|) and taking 
the large time (f3) limit, we obtain the relations between tunneling propaga- 
tors and the level shifts Aei, A €2 respectively, where eo is the corresponding 
ground state energy for both Hq and Hq 2 \ i.e. 

u u ($2+i) e~ 2 ^ sinh(2/3A ei ) 

= e-^- 2 ^ 1 "^^ (<P f = $2+1,/^ = (36) 

and 

^0 ^0 ($2) e- 2/3eo sinh(2/3Ae 2 ) 

= e-'-W**- 1 -*® (0/ = = $2-1,-/3) (37) 

In the following we evaluate the Feynman propagators about the instanton 
trajectories (jyp and (p^ and compare the results with eqs.(^) and ([H]) 
to determine the level shifts Ae l5 At2. The propagators K^\i = 1,2, are 
defined by eqs.(|9|) to (|TT|). We expand 4> about the classical trajectories <f>^\ 

rh = ,A« + v <t) (38) 



c 



where r/^ denotes the small fluctuation about the classical trajectory (jy^ 
with fixed end points such that 

7^(0,) =77« (^ = (39) 

Substituting eq. (|38|) into the propagator eq.(H) and keeping terms up to and 
including those quadratic in 77 W for the one-loop approximation, one has 

K {i) = g-^/W (40) 
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where 



/(*) 



(i)-|_-ASW 



/,,(*) (^)=0 

is the fluctuation functional integral with fluctuation action 



where 



M® = - 
2 



(41) 

(42) 
(43) 



is the operator of the second variation of the action, and 

v-W) - ^ 



We let {V^} be the set of normalized eigenf unctions of the operator MW 
with eigenvalues {-E 1 ^} such that 



We may expand 77 W in terms of 



(0 



Then the functional integral can be formally expressed as 



/(«) 



9c 



(0 



n, 

n=0 •C'i 



7T ,1 
12 



(0- 



(44) 



(45) 



(46) 



The first factor on the right is the Jacobi determinant of the transformation 
fl45|) . Because of the time translation symmetry of the action there exists 
naturally a zero mode of MW (with = 0). This zero mode is seen to be 
just the time derivative of the classical solution, i.e. 



(47) 



dr 



where is a normalization constant. The fluctuation functional is diver- 
gent, as seen from eq.(HH). This well known problem of the quantization 
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about a classical solution can be cured by the Faddeev-Popov (FP) tech- 
nique [14] or in a more systematic way using the BRST transformation [|i~5| . 
The result is that the integration over the zero mode variable, dco, is replaced 
by an integration over the collective coordinate, i.e. the instanton position 



g?To which leads to a factor 2(3 and an FP determinant [15] y D given by 



D 



ffdr = 2K X S® 



One then has 
where 



/« = 20j2K 1 sPlH ) 



r(*) 



dr] 



(0 



(0 



7T ,1 
I 2 



(48) 



(49) 



(50) 



which is not easy to evaluate directly although it is divergence free. We now 
evaluate jW in an alternative way known as the shift method [15 1 in path 
integral theory. Introducing the transformation 



i7«(r)=y«(r) + ^( T ) 



c 1 yW^Orfr' 



(51) 



which converts the functional integral 1^ into a Gaussian integral. We obtain 



J0 



■ 1 ^ 
Z7T 



(*)/ 



r/r 



(52) 



-/3 m($°) 2 _ 

Substituting the explicit formula for 0^ which is the time derivative of 



eqs. (p3|) , (24j), we find the functional fluctuation integrals are the same for 
both instantons ^ and in the large time limit, i.e. (3 — > oo, so that 



= j(2) 



^7T 



(1 - a )*A4s 



i.i i 

2 



(53) 



To obtain the desired result the fluctuation functional integral expression 
may be reexpressed by comparison of eqs. fl49|), fl50|) and ( |53"D as 



J (i) = 2P\J2K 1 S, 



7T 



2tt 



1 — a )4A 4 s 2 



(54) 
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This expression involves explicitly the zero mode eigenvalue Eq which van- 
ishes when (3 — > oo. We can find the explicit time dependence of Eq\/3) 
by using a socalled boundary perturbation method |], which leads to the 
formula 



rn 



-IE, 



{4>f) 2 dr (55) 



where 



\P) 



dr 



HP) = 



d 2 #(r) 



T=0 



dr 2 



T=(3 



After evaluation of both sides of eq.(]55f) we again take the large time limit. 
The time dependence of the zero mode eigenvalue is then found to be 



E 



(i) _ 4A(1-Q! 2 )2g 2 tj ^_ 2VCT 



5, 



(56) 



Substituting this expression into eq. (|54D the desired result of the Feynman 
propagator due to tunneling of one instanton is found to be 



7T 



(57) 



The subscript (1) denotes the one instanton contribution to the propagator. 

The path integral quantization requires a sum over all possible paths. In 
the present case the paths fall into an infinite number of classes consisting 
of one instanton plus an arbitrary number of instanton-anti-instanton pairs 
for a given time interval, i.e. 



(58) 



n=0 



where K$ n+ x\ is the Feynman kernel of the contribution of one instanton 
plus n pairs. The propagator with n pairs is found to be 



JX {2n+l) 



2n+l 



— (1 - a 2 )*s 2 uJo\4 

7T2 



(2n + 1)! 



2n+l 



(59) 
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The final result of the tunneling propagator is 

K®= ^"^V ^^sinh 
V2^ 



71 



(60) 



We have seen that the quantum fluctuation part of the tunneling propagator 
is the same for both types of instantons (f)^ and (j)^>. The only difference 
comes from the classical action. Comparing eq.flSUD with eqs. (|36"D and ([37]) 
we obtain the corresponding level shifts Aei,i = 1,2, 



Aej = P(a, s, A)e 



(61) 



The prefactor is given by 



P(a, s, A) = -^(1 - a 2 ) 4 s 2 cuoA 4 

'71 



(62) 



6. Quenching of macroscopic quantum coherence 



The small ferromagnetic particle described by the Hamiltonian of eq.(|I|) 
in the limit of large spin s has two degenerate easy directions in the absence of 
an external magnetic field, corresponding to the two degenerate ground states 
at 4> = 0, Ti. Quantum tunneling mixes the degenerate ground states and the 
total moment of the small ferromagnetic particle may resonate between the 
degenerate easy directions, a phenomenon known as macroscopic quantum 
coherence (MQC). It has been pointed out that the MQC is quenched for 
half integer spin s, and that the quenching of MQC is seen as a consequence 



of Kramer's degeneracy ||17|| . In earlier work |T(| we demonstrated that the 
quenching of MQC should be understood as the quenching of an energy band, 
and is due to the destructive interference of topological phases of the tunnel- 
ing amplitude in eq.(^). In the present case with a small applied magnetic 
field the easy directions are shifted to = sin -1 a and tt — sin -1 a. The 
quantum tunneling takes place either through the small barrier with angle 
difference tt — 2 sin -1 a or through the larger barrier with angle difference 
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7r + 2 sin^a. The total magnetic moment of the small ferromagnetic par- 
ticle may also resonate between the degenerate easy directions. It may be 
interesting to see the effect of the topological phase on MQC. To this end we 
rewrite the energy spectrum of eq.(|32[) as 



E = e + ( Aei + Ae 2 ) cos [(£ + s)tt] cos 2(£ + s) sin 1 a 
-(Aei - Ae 2 ) sin [(£ + s)%] sin [2(f + s) sin" 1 a 



(63) 



For a weak magnetic field a << 1 . Remembering that £ can assume only 
either of the two values and 1 in the first Brillouin zone, the energy spectrum 
can be rewritten in terms of the explicit level shifts of eqs.(^) and the 
classical actions of eqs. (|26|) and ( p7|) as 

E = e + 2P(a, s, A)e~ 2sv/X j cos[(s + f )7r] cos(2sa) cosh(sv / Aavr) 



+ sin[(s + £)7r] sin(2sa) sinh(sv^Aa7r) 



(64) 



In view of the second term the energy band cannot be quenched for half 
integer spin s, since Kramer's degeneracy is removed by an external magnetic 
field except when sy/Xan is negligibly small. There is, in fact, another type 
of quenching of MQC here for both integer spin s if 



2as 



n + - | 7T 



and half integer spin s if 



2as = mi 



(65) 



(66) 



where n is an integer. This quenching of MQC is irrelevant to Kramer's 
degeneracy. 



7. Conclusions 

In the above we presented a first report on the application of the instan- 
ton method to the analysis of quantum tunneling in a periodic potential 
with asymmetric twin barriers. There are two types of instanton solutions 



15 



corresponding to the two barriers of different sizes. The instanton method 
alone can only yield the separate level shifts due to tunneling. The energy 
spectrum, eq. (|3"2]) , and finally eqs.(^) and here first derived and calcu- 
lated in the one-loop approximation, were obtained with the help of Bloch 
theory and involve both level shifts. A new kind of quenching of MQC is also 
discussed. 
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Figure Caption 

The periodic potential with asymmetric twin barriers and the instanton tra- 
jectories 



18 



